SARS-CoV-2 spike protein diversity at an intra-host level, among SARS-CoV-2 infected individuals in South Africa, 2020 to 2022.

Intra-host diversity studies are used to characterise the mutational heterogeneity of SARS-CoV-2 infections in order to understand the impact of virus-host adaptations. This study investigated the frequency and diversity of the spike (S) protein mutations within SARS-CoV-2 infected South African individuals. The study included SARS-CoV-2 respiratory samples, from individuals of all ages, received at the National Health Laboratory Service at Charlotte Maxeke Johannesburg Academic hospital, Gauteng, South Africa, from June 2020 to May 2022. Single nucleotide polymorphism (SNP) assays and whole genome sequencing were performed on a random selection of SARS-CoV-2 positive samples. The allele frequency (AF) was determined using TaqMan Genotyper software for SNP PCR analysis and galaxy.eu for analysis of FASTQ reads from sequencing. The SNP assays identified 5.3% (50/948) of Delta cases with heterogeneity at delY144 (4%; 2/50), E484Q (6%; 3/50), N501Y (2%; 1/50) and P681H (88%; 44/50), however only heterogeneity for E484Q and delY144 were confirmed by sequencing. From sequencing we identified 9% (210/2381) of cases with Beta, Delta, Omicron BA.1, BA.2.15, and BA.4 lineages that had heterogeneity in the S protein. Heterogeneity was primarily identified at positions 19 (1.4%) with T19IR (AF 0.2-0.7), 371 (92.3%) with S371FP (AF 0.1-1.0), and 484 (1.9%) with E484AK (0.2-0.7), E484AQ (AF 0.4-0.5) and E484KQ (AF 0.1-0.4). Mutations at heterozygous amino acid positions 19, 371 and 484 are known antibody escape mutations, however the impact of the combination of multiple substitutions identified at the same position is unknown. Therefore, we hypothesise that intra-host SARS-CoV-2 quasispecies with heterogeneity in the S protein facilitate competitive advantage of variants that can completely/partially evade host's natural and vaccine-induced immune responses.

Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages.

SARS-CoV-2 lineages and variants of concern (VOC) have gained more efficient transmission and immune evasion properties with time. We describe the circulation of VOCs in South Africa and the potential role of low-frequency lineages on the emergence of future lineages. Whole genome sequencing was performed on SARS-CoV-2 samples from South Africa. Sequences were analysed with Nextstrain pangolin tools and Stanford University Coronavirus Antiviral & Resistance Database. In 2020, 24 lineages were detected, with B.1 (3%; 8/278), B.1.1 (16%; 45/278), B.1.1.348 (3%; 8/278), B.1.1.52 (5%; 13/278), C.1 (13%; 37/278) and C.2 (2%; 6/278) circulating during the first wave. Beta emerged late in 2020, dominating the second wave of infection. B.1 and B.1.1 continued to circulate at low frequencies in 2021 and B.1.1 re-emerged in 2022. Beta was outcompeted by Delta in 2021, which was thereafter outcompeted by Omicron sub-lineages during the 4th and 5th waves in 2022. Several significant mutations identified in VOCs were also detected in low-frequency lineages, including S68F (E protein); I82T (M protein); P13L, R203K and G204R/K (N protein); R126S (ORF3a); P323L (RdRp); and N501Y, E484K, D614G, H655Y and N679K (S protein). Low-frequency variants, together with VOCs circulating, may lead to convergence and the emergence of future lineages that may increase transmissibility, infectivity and escape vaccine-induced or natural host immunity.

Efficacy of Primary Series AZD1222 (ChAdOx1 nCoV-19) Vaccination Against SARS-CoV-2 Variants of Concern: Final Analysis of a Randomized, Placebo-controlled, Phase 1b/2 Study in South African Adults (COV005)

COVID-19 vaccine efficacy (VE) has been observed to vary against antigenically distinct SARS-CoV-2 variants of concern (VoC). Here we report the final analysis of VE and safety from COV005: a phase 1b/2, multicenter, double-blind, randomized, placebo-controlled study of primary series AZD1222 (ChAdOx1 nCoV-19) vaccination in South African adults aged 18–65 years. South Africa's first, second, and third waves of SARS-CoV-2 infections were respectively driven by the ancestral SARS-CoV-2 virus (wild type, WT), and SARS-CoV-2 Beta and Delta VoCs. VE against asymptomatic and symptomatic infection was 90.6% for WT, 6.7% for Beta and 77.1% for Delta. No cases of severe COVID-19 were documented ahead of unblinding. Safety was consistent with the interim analysis, with no new safety concerns identified. Notably, South Africa's Delta wave occurred ≥ 9 months after primary series vaccination, suggesting that primary series AZD1222 vaccination offers a good durability of protection, potentially due to an anamnestic response. Clinical trial identifier: CT.gov NCT04444674

Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing.

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

Molecular Epidemiology and Diversity of SARS-CoV-2 in Ethiopia, 2020-2022.

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499,000 cases, and ~7500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was outcompeted by the Delta variant. Notably, Ethiopia's lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia's role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

Omicron BA.4/BA.5 escape neutralizing immunity elicited by BA.1 infection.

SARS-CoV-2 Omicron (B.1.1.529) BA.4 and BA.5 sub-lineages, first detected in South Africa, have changes relative to Omicron BA.1 including substitutions in the spike receptor binding domain. Here we isolated live BA.4 and BA.5 viruses and measured BA.4/BA.5 neutralization elicited by BA.1 infection either in the absence or presence of previous vaccination as well as from vaccination without BA.1 infection. In BA.1-infected unvaccinated individuals, neutralization relative to BA.1 declines 7.6-fold for BA.4 and 7.5-fold for BA.5. In vaccinated individuals with subsequent BA.1 infection, neutralization relative to BA.1 decreases 3.2-fold for BA.4 and 2.6-fold for BA.5. The fold-drop versus ancestral virus neutralization in this group is 4.0-fold for BA.1, 12.9-fold for BA.4, and 10.3-fold for BA.5. In contrast, BA.4/BA.5 escape is similar to BA.1 in the absence of BA.1 elicited immunity: fold-drop relative to ancestral virus neutralization is 19.8-fold for BA.1, 19.6-fold for BA.4, and 20.9-fold for BA.5. These results show considerable escape of BA.4/BA.5 from BA.1 elicited immunity which is moderated with vaccination and may indicate that BA.4/BA.5 may have the strongest selective advantage in evading neutralization relative to BA.1 in unvaccinated, BA.1 infected individuals.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Identification of SARS-CoV-2 Omicron variant using spike gene target failure and genotyping assays, Gauteng, South Africa, 2021.

The circulation of Omicron BA.1 led to the rapid increase in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases in South Africa in November 2021, which warranted the use of more rapid detection methods. We, therefore, assessed the ability to detect Omicron BA.1 using genotyping assays to identify specific mutations in SARS-CoV-2 positive samples, Gauteng province, South Africa. The TaqPath™ COVID-19 real-time polymerase chain reaction assay was performed on all samples selected to identify spike gene target failure (SGTF). SARS-CoV-2 genotyping assays were used for the detection of del69/70 and K417N mutation. Whole-genome sequencing was performed on a subset of genotyped samples to confirm these findings. Of the positive samples received, 11.0% (175/1589) were randomly selected to assess if SGTF and genotyping assays, that detect del69/70 and K417N mutations, could identify Omicron BA.1. We identified SGTF in 98.9% (173/175) of samples, of which 88.0% (154/175) had both the del69/70 and K417N mutation. The genotyped samples (45.7%; 80/175) that were sequenced confirmed Omicron BA.1 (97.5%; 78/80). Our data show that genotyping for the detection of the del69/70 and K417N coupled with SGTF is efficient to exclude Alpha and Beta variants and rapidly detect Omicron BA.1. However, we still require assays for the detection of unique mutations that will allow for the differentiation between other Omicron sublineages. Therefore, the use of genotyping assays to detect new dominant or emerging lineages of SARS-CoV-2 will be beneficial in limited-resource settings.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Emergence and phenotypic characterization of the global SARS-CoV-2 C.1.2 lineage.

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The initial C.1.2 detection is associated with a high substitution rate, and includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 variants of concern or variants of interest. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta show high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021.

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

T cell responses to SARS-CoV-2 spike cross-recognize Omicron.

The SARS-CoV-2 Omicron variant (B.1.1.529) has multiple spike protein mutations1,2 that contribute to viral escape from antibody neutralization3-6 and reduce vaccine protection from infection7,8. The extent to which other components of the adaptive response such as T cells may still target Omicron and contribute to protection from severe outcomes is unknown. Here we assessed the ability of T cells to react to Omicron spike protein in participants who were vaccinated with Ad26.CoV2.S or BNT162b2, or unvaccinated convalescent COVID-19 patients (n = 70). Between 70% and 80% of the CD4+ and CD8+ T cell response to spike was maintained across study groups. Moreover, the magnitude of Omicron cross-reactive T cells was similar for Beta (B.1.351) and Delta (B.1.617.2) variants, despite Omicron harbouring considerably more mutations. In patients who were hospitalized with Omicron infections (n = 19), there were comparable T cell responses to ancestral spike, nucleocapsid and membrane proteins to those in patients hospitalized in previous waves dominated by the ancestral, Beta or Delta variants (n = 49). Thus, despite extensive mutations and reduced susceptibility to neutralizing antibodies of Omicron, the majority of T cell responses induced by vaccination or infection cross-recognize the variant. It remains to be determined whether well-preserved T cell immunity to Omicron contributes to protection from severe COVID-19 and is linked to early clinical observations from South Africa and elsewhere9-12.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch.

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Sixteen novel lineages of SARS-CoV-2 in South Africa.

The first severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in South Africa was identified on 5 March 2020, and by 26 March the country was in full lockdown (Oxford stringency index of 90)1. Despite the early response, by November 2020, over 785,000 people in South Africa were infected, which accounted for approximately 50% of all known African infections2. In this study, we analyzed 1,365 near whole genomes and report the identification of 16 new lineages of SARS-CoV-2 isolated between 6 March and 26 August 2020. Most of these lineages have unique mutations that have not been identified elsewhere. We also show that three lineages (B.1.1.54, B.1.1.56 and C.1) spread widely in South Africa during the first wave, comprising ~42% of all infections in the country at the time. The newly identified C lineage of SARS-CoV-2, C.1, which has 16 nucleotide mutations as compared with the original Wuhan sequence, including one amino acid change on the spike protein, D614G (ref. 3), was the most geographically widespread lineage in South Africa by the end of August 2020. An early South African-specific lineage, B.1.106, which was identified in April 2020 (ref. 4), became extinct after nosocomial outbreaks were controlled in KwaZulu-Natal Province. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify new lineages and inform measures to control the spread of SARS-CoV-2. Such genomic surveillance presented in this study has been shown to be crucial in the identification of the 501Y.V2 variant in South Africa in December 2020 (ref. 5).

Transmission dynamics of SARS-CoV-2 within-host diversity in two major hospital outbreaks in South Africa.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes acute, highly transmissible respiratory infection in humans and a wide range of animal species. Its rapid global spread has resulted in a major public health emergency, necessitating commensurately rapid research to improve control strategies. In particular, the ability to effectively retrace transmission chains in outbreaks remains a major challenge, partly due to our limited understanding of the virus' underlying evolutionary dynamics within and between hosts. We used high-throughput sequencing whole-genome data coupled with bottleneck analysis to retrace the pathways of viral transmission in two nosocomial outbreaks that were previously characterised by epidemiological and phylogenetic methods. Additionally, we assessed the mutational landscape, selection pressures, and diversity at the within-host level for both outbreaks. Our findings show evidence of within-host selection and transmission of variants between samples. Both bottleneck and diversity analyses highlight within-host and consensus-level variants shared by putative source-recipient pairs in both outbreaks, suggesting that certain within-host variants in these outbreaks may have been transmitted upon infection rather than arising de novo independently within multiple hosts. Overall, our findings demonstrate the utility of combining within-host diversity and bottleneck estimations for elucidating transmission events in SARS-CoV-2 outbreaks, provide insight into the maintenance of viral genetic diversity, provide a list of candidate targets of positive selection for further investigation, and demonstrate that within-host variants can be transferred between patients. Together these results will help in developing strategies to understand the nature of transmission events and curtail the spread of SARS-CoV-2.

Detection of a SARS-CoV-2 variant of concern in South Africa.

Continued uncontrolled transmission of SARS-CoV-2 in many parts of the world is creating conditions for substantial evolutionary changes to the virus1,2. Here we describe a newly arisen lineage of SARS-CoV-2 (designated 501Y.V2; also known as B.1.351 or 20H) that is defined by eight mutations in the spike protein, including three substitutions (K417N, E484K and N501Y) at residues in its receptor-binding domain that may have functional importance3-5. This lineage was identified in South Africa after the first wave of the epidemic in a severely affected metropolitan area (Nelson Mandela Bay) that is located on the coast of the Eastern Cape province. This lineage spread rapidly, and became dominant in Eastern Cape, Western Cape and KwaZulu-Natal provinces within weeks. Although the full import of the mutations is yet to be determined, the genomic data-which show rapid expansion and displacement of other lineages in several regions-suggest that this lineage is associated with a selection advantage that most plausibly results from increased transmissibility or immune escape6-8.